Work exchange between quantum systems: the spin-oscillator model 
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With the development of quantum thermodynamics it has been shown that relaxation to thermal 
equilibrium and with it the concept of heat flux may emerge directly from quantum mechanics. This 
happens for a large class of quantum systems if embedded into another quantum environment. In 
this paper, we discuss the complementary question of the emergence of work flux from quantum 
mechanics. We introduce and discuss two different methods to assess the work source quality of a 
system, one based on the generalized factorization approximation, the other based on generalized 
definitions of work and heat. By means of those methods, we show that small quantum systems can, 
indeed, act as work reservoirs. We illustrate this behavior for a simple system consisting of a spin 
coupled to an oscillator and investigate the effects of two different interactions on the work source 
quality. One case will be shown to allow for a work source functionality of arbitrarily high quality. 



I. INTRODUCTION 

Thermodynamics is a theory of impressive success and 
a wide range of applicability. This is true despite its ori- 
gin as a purely phenomenological theory in the late 16th 
and 17th century. It took about two hundred years until 
Boltzmann and Maxwell gave a foundation of thermody- 
namics in terms of statistical mechanics based on classical 
mechanics. It was only in recent years that a new ap- 
proach to the foundation of thermodynamics from a mi- 
croscopic theory has received increased attention: quan- 
tum thermodynamics, the derivation of thermodynamics 
from quantum theory. In the seminal work of Gemmer 
et al. pi H, E| , the emergence of thermodynamic behav- 
ior from quantum mechanics for a wide class of quantum 
systems has been established. Other recent papers have 
extended and clarified the results of the aforementioned 
authors 0,0,0. 

The interest in a quantum approach to thermodynam- 
ics is two-fold: First, it promises a deeper understand- 
ing of thermodynamic core concepts like relaxation, irre- 
versibility, heat, work, and might elucidate connections 
between those and concepts known from quantum me- 
chanics, e.g., entanglement. Second, such an approach 
should help to find generalizations for the mentioned con- 
cepts for non-equilibrium, finite systems and strong in- 
teraction. 

The definition of work in quantum systems has been 
discussed in various papers 0,0, Q and have been app lied 
successfully to quantum heat engines [Tol [TTl Il3l [T3| . 
Still, all those investigations typically deal with quan- 
tum systems that are subject to driving by means of a 
time-dependent Hamilton operator of the system. Thus, 
the identification and definition of work is determined a 
priori by relating it to the presence of classical driving, 
while no microscopic derivation of the concept is given. 
In the present paper we deal with closed, finite quantum 
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systems. For such systems the functionality of cooling 
has already been investigated (see, e.g. p~5[). Here we 
are interested in the question under what conditions a 
quantum system coupled to another can exert the effect 
of a classical driver over the other system and thus be 
identified with a reversible work source. By our new ap- 
proach based on a complete quantum modelling of the 
work source, we are able to show that classical driving 
and therefore work is not a concept bound to macroscopic 
devices. 

The paper is organized as follows. In Sec. [Til we present 
the factorization approximation (FA) and its generaliza- 
tion to the case of semi-mixed factorizing initial states. 
It is shown that the applicability of the FA allows one to 
identify quantum systems as classical drivers. In the sub- 
sequent Sec. lIVl we deal with the question of work source 
quality definition. Based on the previous section, we in- 
troduce a measure inspired by the FA and establish its 
connection to work source functionality. In addition, we 
develop another work reservoir quality measure based on 
considerations of work and heat fluxes, if an appropriate 
definition of those is given. The definition of work and 
heat flux we chose to use throughout the paper is taken 
from [l6| and outlined in Sec. IIIII In Sec. |V] we present 
the spin-oscillator model and its properties for two types 
of interactions. This simple quantum model is then used 
to illustrate the implementation of work sources of ar- 
bitrarily high quality in quantum mechanics in Sec. I VII 
and to discuss the limits of two different work function- 
ality quality measures in Sec. IVIII where we also give an 
idea of the overall work source behavior of the second 
type of the model. Finally, we summarize our results in 
Sec. EM 



II. EFFECTIVE DYNAMICS: 
TIME-DEPENDENT DRIVING 

The factorization approximation (FA) has been thor- 
oughly discussed, e.g., in p~3| and p~8[. We will therefore 
only summarize the basic statements and give a general- 
ization to the result of 1171] - 
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In its form stated in [17(, the FA reads as follows. Let 
us consider a bipartite quantum system with Hamilton 
operator 



H = Hi + H i 



■Ho 



(1) 



acting on the joint Hilbert space H.12 — Tii <8> "H%. The 
operators H\ and Hi act on the respective local Hilbert 
spaces 7ii and 7i 2 only. Let the initial state factorize, 



i.e. 



|*(0)) = |*i(0))® 1*2(0)). 



(2) 



After some time t, the total state of the system either 
given by |*(i)) or its density matrix pit) = |*(t))(*(t)| 
gives rise to the reduced states of the two subsystems 
fk(t) = tr 2 (|*(t))<M>(t)|) and p 2 (t) = tri(]*(*)><*(*)|). 
Although the state was assumed to factorize initially, 
in general, the subsystem states no longer will be pure 
states due to entanglement introduced by the interac- 
tion between the subsystems. As the total state is 
pure, the subsystem purities P[pi{t)\ — tr{[/5i(i)] 2 } and 
P\pi{t)] — tr{[p2(i)] 2 } are equal at any instant t. Now, 
as long as the purity of the subsystems is close to 1, 
it can be shown that the dynamics of the system are, in 
good approximation, given by the reduced density matri- 
ces p t = \^i(t))(^i(t)\, i = 1,2, where the |*i(t)) obey 
the coupled differential equations 



ifc|*l(t)> = (ffi + (tf 2 (*)l#12|*2(t)>J|*l(t)> (3) 
ift|* 2 (t)) = (#2 + <*l(t)|#12|*l(t)>) |* 2 (*)> (4) 



up to an irrelevant relative phase (for a detailed deriva- 
tion see [I2I )• Obviously, the reduced states of the sys- 
tem evolve under the action of time-dependent effective 
Hamiltonians H\+ H^ s (t) and H 2 + H^ {t), respectively, 
with 



Hf{t) = tr 2 {iT 12 [i®p 2 (t)]} 
Hf{t) = trx{H 17 .[pi{t) ® 1]} 



(5) 
(6) 



and in the present case fij(t) = \^j(t)}{^j(t)\, j = 1,2. 

The above statement can be generalized for the case 
of factorizing semi-mixed initial states, that is, states of 
the form 



p(0)=pi(0)®|* a (0))<* 2 (0)| 



(7) 



If the purity P[p 2 (t)] of the initially pure system 2 re- 
mains close to unity, the dynamics of the system can be 
given approximately by the coupled differential equations 

i%h{t) = [Hi+Hf (t),M*)] (8) 
ifi|* 2 (t)) = (# 2 + £§*(<)) |* a (t)) (9) 

with p 2 (t) = |* 2 (t))(* 2 (t)|. Again, the effect of the sub- 
systems on each other is to induce a time-dependent ef- 
fective Hamiltonian that governs the time evolution of 
the subsystems. For the derivation, see App. [A] 



Here, we would like to stress the fact that, as long as 
the prerequisites for the FA are met, Eqs. ([H [9|) present 
an alternative description of the system dynamics: The 
subsystems can be considered classical drivers for each 
other. It is remarkable that this feature is reciprocal 
and based on (approximate) constancy of the subsystem 
purities (entropies). 

The energy exchanged this way can aptly be called 
"work": Classically one would define the work W im- 
parted over time ts on a Hamiltonian system H (A) with 
A denoting the time-dependent control parameter as fl~oj ] 



/ , dXdH.^. 



(10) 



where z{t) denotes the system's state trajectory in phase 
space. One notes, however, that the energy exchange 
will, in general, be contaminated by contributions violat- 
ing the constancy of local purity. This contamination is 
a characteristic feature of the underlying total (unitary) 
dynamics. Close to thermal equilibrium such a contribu- 
tion would be called heat, &~Q: Work and heat in o pen 
quantum systems are usually defined as 0, H, [TFj, [Tll.lT2T| 



dU = d(H) 



tr(pdH) + ti(Hdp) 

dW SQ 



(11) 



again recognizing the energy exchange in the FA scenario 
as work. 

We emphasize here, that explicitly time-dependent 
Hamiltonians are not part of the fundamental descrip- 
tion of nature as given by quantum mechanics. There- 
fore, there is no way how they could come about save by 
an effective description of a system like the FA. If one 
denied any physical significance of such an effective de- 
scription and hence considered it only a mathematical 
simplification without physical meaning, one obviously 
would have to deny the physical existence of classical 
drivers altogether. This is not a reasonable option. 



III. LEMBAS PRINCIPLE 

The effective dynamics according to Eqs. (|8l9p allows 
for an intuitive approach to the concept of work: In gen- 
eral, however, only approximately; the deviations remain 
unquantified. 

Here, the LEMBAS approach [TJI comes into play 
based on the following ideas: First, choose a partitioning 
of the total isolated system into system of interest (1) 
and its environment (2) and consider the exact local dy- 
namics of the system (1). The state of the total system 
is 



p = pi ® p 2 + C\ 



(12) 



where pj are the respective reduced density operators. 
Then, the exact (effective) Liouville-von Neumann equa- 
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tion for subsystem (1) can be written as 



H 1 +Hf{t),p 1 (jt)\+£f\p{t)] (13) 



cfw)] 



with the superoperator 

-ih-hv 2 {[H 12 ,c 12 }}. 

The local energy H[ is defined now based on consider- 
ations how the local system would appear to an experi- 
menter ("local effective measurement basis", LEMBAS). 
There is some ambiguity in the procedure, but it has 
proven useful in [l6| to choose 



(14) 



where Hf(t) = Hf a (t) + Hf b (t) and Hfjt) is the part 

of Hf s (t) that commutes with H\. 

The final step is to discriminate energy changes of the 
system based on whether they change the local von Neu- 
mann entropy 5*i or not, that is whether they are of co- 
herent (work) or incoherent origin (heat). This leads to 
the following formulas for heat- and work-flux for any 
partitioning and any H[: 

W(t) = tr{A 1 cff a (t)p 1 (i)- i [^W,^i cff 6 W]PiW}(15) 

Q(t) = tv{H{(t)£f[p(t)]}- (16) 

How do these generalized definitions connect to their 
thermodynamic analogues? In the thermodynamic limit, 
that is, close to the thermodynamic equilibrium, for in- 
finitely sized systems and weak couplings, the von Neu- 
mann entropy of the respective subsystem and its ther- 
modynamic entropy coincide and the LEMBAS defini- 
tions of work and heat blend in with their thermodynamic 
counterparts. 

But also in far from equilibrium situations, the LEM- 
BAS definitions can be associated with work and heat in 
the following sense: We know from the results of quantum 
thermodynamics [3, 0, Q that thermodynamic behavior 
of a system can be seen to result from an embedding in 
an environment, which by itself needs not to be and usu- 
ally is not thermodynamic (in equilibrium, infinite, weak 
coupling). Thus, validity of thermodynamic concepts is 
not a property of the total system but has to do with 
whether or not the system of interest is influenced by its 
environment in such a way that thermodynamic proper- 
ties emerge, which is a purely local consideration. The 
LEMBAS definitions take this concept to the extreme in 
the sense that they state that "what locally has a work 
effect Hf s (t), is work" and "what locally has a heat effect 
Cf[p ( t)}, is heat" even for non-thermodynamic (in the 
classical sense), far from equilibrium situations. Mak- 
ing the distinction in this way is justified by the fact 
that classical driving can be unambiguously identified as 
work even in the thermodynamic sense and, therefore, 
any effect £° ff [/>(£)] not related to work is identified as 
heat. 

Finally, we note that the LEMBAS definitions retain 
the properties that 



1. work is energy exchange due to changing param- 
eters of the Hamilton operator that describes the 
system; 

2. heat is energy exchange associated with change of 
entropy, although here a generalized definition of 
entropy is to be used. 



IV. MEASURES OF WORK SOURCE QUALITY 
A. Work reservoir 

An ideal work reservoir can be defined as a system ex- 
changing energy only in the form of work. It is obvious 
that this definition is too restrictive for the classification 
of realistic models, that is, models involving finite size, 
finite interaction and limited control. No realistic model 
can comply to the idea of such an ideal work source as 
even arbitrary small but finite deviations from this ide- 
alized concept would lead to a rejection of a model as a 
work source. Additional complications arise due to the 
fact that we have to consider processes, the properties of 
which may change with time. 

Thus, there is need for a more differentiated measure of 
work reservoir functionality. In a non ideal world, special 
attention is to be paid to the definition and quantification 
of the quality of a work reservoir to be able to compare 
and to draw conclusions on justified grounds. 

Basically, one can distinguish two types of measures 
depending on whether they refer to a single point in time 
or to a (finite or infinitely large) interval of time. Wc like 
to refer to them as instantaneous and integral measures 
and our main interest lies on the integral ones, defined 
with respect to some finite time interval (again because 
under realistic condition it is not expected that a system 
can be a work source for all times) . 

In the following section, we present two different ap- 
proaches to the problem based on two distinct physical 
reasonings. 

B. Purity based measure 

Comparing Eq. (TL5|) to Eq. (jHJl, one realizes that the 
applicability of the FA is equivalent to a vanishing £° ff . 
Thus, if the total system was initially in a semi-mixed 
state, £f ff is negligible if P[/5 2 (i)] « 1. In this sense, 
P\p2{t)] is a measure of work reservoir functionality. The 
closer it is to 1, the smaller £\ s has to be and the less 
energy may be exchanged as heat instead of work. Note, 
that acting as a work reservoir is a reciprocal property, 
i.e., each subsystem acts on its partner in an analogous 
way. This is in perfect agreement with what we know 
from thermodynamics. If we have two systems under- 
going a process during which only work is exchanged 
between them, both systems obviously act as work reser- 
voirs for each other although we may imagine one system 
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to be the gas filling a box and the other system to be the 
piston capping the box and being connected to a spring. 

At first glance, the purity therefore seems to be a good 
candidate for assessing work source quality: It is an easy 
quantity to compute - even analytically - and by its con- 
nection to the FA, the physical reasoning is clear. 

However, as clear as the ideal situation with P[pa(t)] = 
1 is, it is unclear to give a quantitative interpretation 
for purities lower than unity because there is neither an 
obvious relation between P[p2(t)] and £f ff nor between 
£f and the quality of the work reservoir functionality. 
Moreover, it is expected that the same purity decrease 
for different systems, especially of different size, has to 
be weighted differently. Thus, any concrete choice of a 
minimum purity beyond which a system will be accepted 
as a work reservoir will remain somewhat arbitrary and 
difficult to compare with other systems' purity behavior. 
If such a purity threshold was given, the respective sys- 
tem could be considered as a work source for any time 
interval during which the purity stays above the given 
threshold. 

As will become evident in Sec. IV Bl there is another 
problem besides the arbitrary definition of the threshold 
when using this measure: The decrease in purity is linked 
to the size of C\ only. Thus, the purity does not contain 
any information about the relative effects of Hf s (t) and 
£° ff . Since the former is related to work and the latter 
to heat, a comparison of both in terms of their effect on 
the energy of the system is in general expected to be an 
important part of the assessment of work source quality. 



Work and heat flux based measure 



integrate and then normalize the result defining 



We introduce the ratio 
r(i) 



\W(t)\ 



(17) 



\W(t)\ + \Q(t)\ 
which has the following convenient properties: 

• r (t) = 1 W(t) =/= A Q(t) = 0: ideal work source 

• r(t) = & W(t) = A Q(t) ^ 0: ideal heat source 

Provided there is energy exchange at all (i.e. not both, 
W, Q are zero), r is well behaved. As we took separate 
moduli in the denominator, there can be no compensa- 
tion due to opposite sign. 

Based on this instantaneous measure, we can develop 
an integral measure for finite time intervals [io>ii]- Di- 
rectly integrating over r(t) is not an option for this would 
completely ignore the time-dependence of the total of the 
absolute fluxes and therefore the necessary weighting of 
r. It is straightforward to apply the necessary weight, 



R(ti,t ) := 



Jr(t)(\W(t)\ + \Q(t)\ 



dt 



j(\W(t)\ + \Q(t)\ 

to V 

J\W(t)\dt 



dt 



J(\W(t)\ + \Q(t)\ 



(18) 



dt 



Defining the quantities 



W(h,t ) 



\W(t)\dt, Q(h,t ) 



\Q(t)\dt (19) 



in 



we can rewrite Eq. (fT8"|) in the form of Eq. (jTTJ) as 



R(tx,to) 



W(ti,to) 



W(h,t ) + Q(h,t ) 



(20) 



This integral measure has the same special points like the 
instantaneous measure with the following interpretations: 

• R(t) = 1 & Q(t) = for all t € [i ,ii] and W(t) ^ 
for some t G [io,ii] : ideal work source 

• R(t) = & W(t) = for all t € [t , h] and Q(t) + 
for some t £ [^oi^i] : ideal heat source 

We stress the fact here that a measure based on the in- 
tegrated work W(ti,t ) = f t * W(t)dt and the (analo- 
gously) integrated heat is not able to accomplish such 
precise assessment of the work source quality: For oscil- 
lating fluxes, e.g., Q(ti,t ) might reach for some inter- 
val, although during the time interval there might have 
flown vast amounts of heat. By employing the integrals 
of the absolute fluxes in the chosen definition, we achieve 
a much stronger statement about the quality of a system. 

Finally, let us note that there is also a drawback to this 
measure, namely the difficulty of calculating it because 
of the integration over the absolute values of the fluxes. 



V. APPLICATION: 



SPIN-OSCILLATOR MODEL 
(SOM) 



We turn now to the description of the model we will use 
to demonstrate the existence of small quantum systems 
that do act as work sources. We illustrate the features 
of FA and the various work measures we have discussed 
above and discuss the model and its properties with spe- 
cial focus on the dynamics of the purity. 

The model is a single spin interacting with a har- 
monic oscillator (spin-oscillator model, SOM). On the 
one hand, the SOM serves as an allusion to a classical 
steam engine with a gas of some temperature (spin) and 



a piston periodically compressing and expanding the gas 
(oscillator). On the other hand, the SOM has been used 
in previous related works as a central element of quantum 
thermodynamic machines .11, 12, 20]. Also, the simplic- 
ity and therefore partially possible analytical treatment 
of the model has further motivated the choice. 
The SOM is defined by the Hamiltonian 



H — -^®z + -Hint 



aa+ 2 



(21) 



where we have set H = 1. We denote the spin and oscilla- 
tor local Hamilton operators as H s and H Q , respectively. 
The eigenstates of H s are |0) and |1) with the respective 
eigenvalues =Fw s /2. The eigenstates of H a are defined as 

{|fc)} with eigenvalues cj (k + 1/2), where k = 0, 1,2, 

We will discuss the z- and the ^-interaction, H lnt = 
{H Z ,H XZ }, where 



Hint = H z = \a z x 

H m t = H xz = X(a z + Ka x )x. 



(22) 
(23) 



For the initial state of the total system, we assume 
that the spin has interacted in the past with some heat 
bath in order to establish a thermal state but now is 
decoupled from said bath (or the bath coupling is so weak 
that its influence may be neglected during the period of 
evolution one is interested in) . The oscillator is prepared 
in a coherent state \a). Thus, the initial state is given as 



P(0) 



c 
1- 



® \a)(a\ 



(24) 



where the spin's state is given in its energy eigenbasis. 
The self-generated process imposed on the spin via cou- 
pling to the oscillator might thus be called "adiabatic"; 
however, due to quantum mechanical interactions the lo- 
cal entropy (purity) will, in general, not be constant, see 
below. 



A. z-interaction (z-SOM) 

Representing the Hamiltonian (|2ip in the eigenbasis of 
the spin, one finds that it has a block-diagonal structure: 



H = 



H — Xx — ~2~ 

H_ 

H 4 



H + Xx + ^f- 



(25) 



The same obviously holds for the time-evolution operator 
and - by the block-diagonal structure of the initial state 
- also for the propagated state p(t) of the total system: 



c|a_(t))(a_(t)| 



(l-c)\a + (t))(a + (t)\ 



(26) 
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FIG. 1: Purity dynamics of the oscillator in the z-SOM for 
the special parameters A = 0.1, c = 0.7, a = 0, m = u s = 

LOo = 1. 



Here we have used the definitions \a±(t)) := U±(t,0)\a) 
and U±(ti,to) — exp [— iH±(ti — to)]- Note that the dy- 
namics of the system are periodic, because both Hamil- 
ton operators H± describe (displaced) harmonic oscil- 
lators with the same frequency lu . Thus, we have 
U (ti+2irmu!~ 1 , t$ + 2itnu)~ 1 ) = U(ti,to) for integer num- 
bers n, m. 

Because of the simple structure of the time-evolution of 
the system, the purity of the oscillator can be computed 
analytically and turns out to be given by 

P[p (t)} = c 2 + (1 - c) 2 + 2c(l - c)|(a_(i)|a + (<))| 2 (27) 



with the time-dependent part 



\(a-(t)\a+(t))\ 2 =exp 



■ sin 



-UJ t 



(28) 



(for the derivation, see Appendix [Bl m is the oscillator 
mass). For pure initial spin states (c = 0,1) we have 
P[p (t)] — 1. An example for P[p (t)] for a mixed initial 
spin state is given in Fig. [T] 

It is easy to see from Eq. (|27I28[) that the minimum 
purity with respect to t and c is 



P" 



where 



1 



[l + exp(-8£)j 



A 2 



Therefore, one has to choose £ — + and thus 
1 - P' nin « 1 



(29) 



(30) 



(31) 



in order to apply the FA. 

We can distinguish two different ways to enforce the 
limit £ -> 0: 

to — * oo, lu = const., A = const. (32) 
uo — ► oo, to = const., \ 2 /lu — const. (33) 
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Their relevance will become clear in Sec. I VI I If one ac- 
cepts the resulting finite P™ ln for some finite £, the local 
coherence-time may be called infinite. 



B. ^-interaction (iz-SOM) 

We discuss now the more complicated case of an inter- 
action of form Eq. (|23|) . This interaction is motivated by 
the following considerations. First, the above case is very 
special in that the minimum purity reached can be con- 
trolled completely by the system parameters. The more 
general xz-interaction case will show that this property is 
dependent on the interaction. Second, the chosen inter- 
action allows us to study the effect of imperfect control 
over the exact form of the interaction as it might be the 
case for a more realistic experimental situation. Finally, 
the xz-interaction exhibits a remarkable diversity of dy- 
namics which serves to illustrate the pros and cons of the 
proposed work reservoir quality measures as well as the 
possibility to realize quantum work sources within the 
given model. 

First, we show that a Hamiltonian 



H = — a z + Asx 



a'a + 



1 



(34) 



with an arbitrary operator s acting on the spin's Hilbert 
space is equivalent to the xz-SOM. This is seen from the 
expansion of s in the operator basis {Qk\k = 0, 1, 2, 3} = 
-^{l,a x ,(Ty,(T z }, which reads 



fe=0 



tr(Qls)Q* 



(35) 



(for details refer to [2J|, pp. 34-49). The 1 term is lo- 
cal to the oscillator and can be absorbed in H Q , while 
the &k terms can always be transformed to the form 
s = X(a z + Ka x ) by the local transformation exp(— i(f>a z ) 
with appropriately chosen real <j). This corresponds to 
a rotation around the z-axis of the spin. The xz-SOM 
discussed here is therefore representative for the whole 
class of Hamiltonians of the form Eq. (fM)) . 



C. Dynamics 

Now, we want to look into the behavior of the system 
for | ac | < 1. To get insight into the dynamics, we invoke a 
rotating wave approximation (RWA). For that purpose, 
we first write the xz-SOM interaction Hamiltonian H xz 
in the interaction picture 

/.n-A\ . „ t / n- a 

H xz oc CT z aexp \i — - — t I + a z a' exp \—i — - — t 

+ A(7_|_a exp(— i At) + A<7_aexp(if2i) 

+ A<r+<r exp(— iflt) + \&-d) exp(iAi)], (36) 



where we have defined :— lu s + lj and A := u> s — uj q . 
By restricting ourselves to the resonant case A = and 
omitting all terms rotating with frequencies and f2/2, 
the xz-SOM Hamiltonian in RWA turns out to be 



H 



RWA 



:&z + g{o+a + <r-a T ) + uj a'a + 



1 



(37) 



in the Schrodinger picture. This is just the Hamiltonian 
i?jc of the Jaynes-Cummings model ( JCM) [22j, |23( with 



Ak 



fJ 



(38) 



and ui = u s = w Q . According to [2J, [23, the RWA is 
accurate as long as g/il <C 1. This condition is met for 
all relevant cases, since we consider in the following a 
situation where the parameters m,uj, and A have been 
chosen such that for n — (z-SOM) the FA holds for all 
times. 

Now let us turn to the interpretation of the result. 
First, we note that by performing the RWA, in particular, 
the <7 Z term of the interaction Hamiltonian is removed. 
Therefore, the xz-SOM in RWA captures the effect of 
the a x interaction alone and, in turn, this means that 
the main dynamics are governed by the a x part of the 
interaction alone. 

Second, the dynamics of the JCM (and therefore of the 
xz-SOM in RWA) scale in time with g _1 . This is most 
clear from the exact time-evolution operator for the JCM 
P- 205) in the case of exact resonance 



cos(gtB) -id) sm{gtA)A^ 1 
-i sm(gtA)A~ 1 a cos(gtA) 



(39) 



with A — V at a + 1 and B = %/ a) a. 

The numerical results of the dynamics of the purity of 
the oscillator are given in Fig. [2] for the case of a coherent 
initial state with one photon in the cavity on average 
(a = 1). The deviations of the numerical results for the 
xz-SOM with and without RWA for three different orders 
of magnitude of g are given in Fig. [3] 

From Fig. [3] one sees that the RWA yields good results 
(less then 10% relative deviation) up to n w 10. This 
shows again that the RWA gives an accurate description 
of the xz-SOM dynamics in agreement with the expecta- 
tion given above. 

It is obvious from Figs. 11121 that the purity behavior of 
the xz-SOM is fundamentally different from the z-case. 
The decrease of the purity due to the additional a x inter- 
action term is several orders of magnitude stronger than 
what is expected from the a z term alone and due to the 
approximate scaling behavior of the xz-SOM, the min- 
imum does not depend on k, as long as k is not zero. 
We conclude from that, that in the presence of an arbi- 
trarily small but non-vanishing a x term the FA and with 
it the work reservoir quality of the oscillator will break 
down in finite time. Reduction of n can only delay the 
breakdown and, thus, if A, k ^ no choice of the other 
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FIG. 2: Purity of the oscillator for numerically exact dynam- 
ics of the xz-SOM in RWA for arbitrary g. The parameters 
are: A = 0.01, a = cj = ui a = 1, c = 0.7. 




FIG. 3: Common logarithm of the absolute deviation 5 of the 
(numerically exact) purity dynamics with and without RWA 
for (top to bottom) g — —0.01, —0.1, —1 corresponding to k = 
\/2, 10\/2, 100\/2 according to Eq. (|38[) with the parameters 



A 



-0.01, m = l,a 



1, c = 0.7. 



model parameters can prevent the breakdown. This re- 
sult is in agreement with the finding in [171 ] for the that 
the coherence time for the JCM depends on the interac- 
tion strength such that weaker interaction leads to longer 
coherence time. 

In the above sense, the work reservoir functionality 
in the given quantum scenario is quite sensitive to the 
quality of control of the interaction between spin and 
oscillator. 



VI. WORK SOURCE QUALITY IN THE z-SOM 

Considering the first case of z-SOM, let us assume that 
one has chosen the parameters of the system such that 
Eq. (|3"Tj) is fulfilled. We can then apply the FA not only 
to describe the dynamics of the system up to any desired 



accuracy but moreover, we get a new level of description 
combined with new physical insight in the properties and 
characteristics of the system. This will be outlined below. 

Applying the FA to the SOM, we find according to 
Eqs. JS]) and ([9]) the following effective coupled equations 



ip s = 
*!*(*)) = 



— + A (£)(*)) cr z ,p s 



H + \(l-2c)x) |*(t)), 



(40) 
(41) 



where we have defined (x)(t) := (^f(t)\x\^(t)). Hence, 
the spin is driven by the oscillator displacement which 
acts like an additional time-dependent magnetic field, 
modulating the spin's Zeeman splitting. On the other 
hand, the spin dynamics lead to a constant displace- 
ment of the oscillator potential. In this sense there is 
asymmetry between the two subsystems: The effective 
Hamiltonian for the oscillator is modified but not time- 
dependent. 

Clearly, this result is in agreement with the result 
which one would obtain from applying the integral mea- 
sure based on W and Q: The latter is zero, since accord- 
ing to Sec. IV Al and Eq. (|40|) . the spin's state does not 
change and the oscillator only exerts classical driving on 
the spin. The local effective energy change of the spin 
resulting from the driving is 100% work. Hence, the os- 
cillator acts as an ideal work source at least in the limits 
discussed in Sec. fVAl Eqs. (|32I33|) . 

However, it has to be noted that the peak-to-peak 
amplitude Aw| ff = A((a;) max — (x) min ) of the effective 
spin splitting is dependent on the system parameters (see 
App.0: 



Alo' 



eff 



2A 



mu! 



-|a + 7|, 



(42) 



where 7 = ^ X 2 / (2muj3) (1 - 2c) = *Jtj2(l - 2c) [cf. 
Eq. f30])]. For the first limit proposed in Eq. l[3"2]). we 
therefore find that both 7^0 and Awf — ► 0, if all 
parameters besides m are kept constant. Thus, the work 
effect induced by the oscillator diminishes more and more 
for increasingly better fulfilled FA. This can be avoided, 
though, by additionally imposing a — > 00, such that 
\a\ 2 /m remains constant, which then defines the split- 
ting's amplitude. This is a classical limit in that the 
mass and average excitation number of the oscillator go 
to infinity. 

There is also a true quantum limit, though, which is 
found to be realized exploiting the second limit given in 
Eq. (J33J). Here, by letting uj — > 00, we enforce that 
£ — > so that the factorization approximation becomes 
exact. However, by requiring X 2 /uj = const., the pref- 
actor of Aw° ff in Eq. (|32"j). 2y / 2\ 2 /(mui ), becomes con- 
stant. Thus, although 7^0, Awf ff retains a finite value 



eff 



2A 



muj 



\a\ = const. 



(43) 



for arbitrary (small but finite) m and a in the limit of 
exact FA. 
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By the preceding reasoning, we conclude that the os- 
cillator is, indeed, a work reservoir for the spin, period- 
ically changing the spin splitting and therefore transfer- 
ring work to/from it. What is special about that finding 
is the fact that a true quantum system (the oscillator 
in the quantum limit of the z-SOM) can be set up as 
an ideal work source and, thus, the work concept is not 
tied to classical devices. Moreover, as long as we fulfill 
Eq. (|3Tj) the oscillator behaves as a work reservoir for any 
time-period. 



VII. WORK SOURCE QUALITY IN THE 

iz-SOM 



In this section, we present and discuss our results for 
the work reservoir behavior in the more general xz-SOM 
presented in Sec. IV Bl with focus on the suitability of the 
work source quality measures proposed in Sec. HVl 

The numerical results used herein have been produced 
with the Mathematica package using the following tech- 
niques: We have computed the time evolution of the 
system by direct diagonalization of the Hamilton oper- 
ator (with a cut-off chosen such that only states with 
occupation probability higher than 10 -6 are included). 
Integration of quantities - where necessary - has been 
performed using the rectangle rule and the error of inte- 
gration has been controlled by crosschecking with results 
for the trapezoidal rule and/or for smaller time steps. 



A. Purity based approach 

We need now to define a lower bound for the purity 
of the oscillator. The a z coupling alone already leading 
to some limited purity loss in the oscillator can be con- 
sidered as sort of a "natural" purity drop, which has to 
be accepted for any system that interacts at all and that 
is present even if the FA is good and the work source 
quality high. 

The work source functionality is considered to fail 
when the purity decrease of the RWA dynamics (caused 
by the a x interaction term) reaches the maximum pu- 
rity drop P™ ln of the a z dynamics alone found for the 
given system parameters (A, oj s , ujq, m). This allows us 
to define a breakdown time t* by P[p (t*)] = P™ m . 

The close connection of the xz-SOM dynamics to the 
JCM dynamics seems to suggest an analytical approach 
based on the standard approximations made to solve the 
JCM (see, e.g.. 



271 ] . Ch. 6): approximation of the occu- 



pation probability of the coherent state with a Gaussian 
and linearization of the spectrum of the JCM around its 
peak. With those approximations a fairly accurate de- 
scription of the JCM's typical collapse and revival be- 
havior of the spin polarization for high initial photon 
numbers (|a| ^> 1) is possible. After the initial collapse 
of polarization, the spin reaches its minimum purity [231 ] 



and the oscillator purity will as well have dropped signif- 
icantly. 

Unfortunately, even in the high photon number limit 
the accuracy of this approach in the relevant time interval 
up to this point of the evolution is insufficient: Typical 
values of the purity drop due to the a z interaction are 
of the order of 10~ 2 , while the error of the mentioned 
approximations is of around the same order during the 
collapse. This renders the application of those approxi- 
mations futile and since a full analytical analysis is much 
too involved, we will only exemplify some results based 
on numerics. 

To this end, we choose the following parameters for 
the xz-SOM: ui — uj s = 1 (resonant case), m = 1, 
A = n = 0.1. In the following, we consider the results of 
two special cases: 

(a) a = 0, c = 0.5 

(b) a = 2, c = 1 

These two examples are drawn from a set of results for 
initial states with parameters a € [0,4] and c £ [0.5, 1.0] 
and have been chosen for they represent in some sense 
extremal cases, that will be seen to illustrate the features 
of the different work source quality measures. A short 
overview about the more general behavior of the xz-SOM 
is given in Sec. IVII CI 

The purity behavior of the examples is shown in Fig. 2] 
The time after which the FA is estimated to fail is roughly 

ft , w 28 and tl,, w 73. Although this means that the 

(a) (&) t. 

second case is expected to exhibit work reservoir func- 
tionality about three times longer than the first case, one 
would conclude from the curves that for both cases, the 
oscillator's work source functionality degrades quickly af- 
ter the initial high purity phase and is virtually absent 
at least for t > 100. 



B. Work/heat based approach 

However, taking a look at the result for the integral 
quality measure R shown in Fig. O one comes to a com- 
pletely different conclusion: In case (a) the oscillator 
starts as a perfect heat source rather than a perfect work 
source and only in the course of time a work source effect 
arises, whereas in case (b) the oscillator is recognized as 
a nearly ideal work source during the whole interval. It 
is astonishing to see that the purity based measure gives 
such a different picture since the reasoning based on the 
FA is valid: During the initial phase, C° is close to 0. 

The reason for this seemingly contradictory character- 
ization becomes evident when examining the results for 
the integrated work and heat, W(t) and Q(t) (see Fig. [5]). 
In both presented cases, the total heat flow in the begin- 
ning of the dynamics is small as expected from the FA 
argument. Also, the heat becomes significant not before 
tT s and t* b y respectively, which again demonstrates the 
strong connection between heat and purity. 
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(a) a = 0, c = 0.5 



200 




200 



FIG. 4: Oscillator purity behavior of two special cases of 
xz-SOM and comparison with minimum purity of z-SOM 
for the given system parameters: Numerical exact result 
(solid line), numerical result with RWA (dashed), minimum 
z-SOM purity (dotted), which - according to Eq. (I29|l - is 
[1 + exp(— 2/25)]/2 w 0.962. The insets show the crossings of 
RWA-purity with minimum z-SOM purity. 
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FIG. 5: Results of the integral work reservoir quality measure 
R(t, 0) for the examples as of Fig. [4] (a) a = 0, c = 0.5 (solid 
line), (b) a — 2,c = 1 (dashed) 
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FIG. 6: Integrated work W{t,0) (solid line) and heat Q{t,0) 
(dashed) for the two chosen examples, as of Fig. [S] The inset 
of case (a) shows Wit, 0) alone. 



However, the work exhibits a completely different be- 
havior for the two cases: In the first case, the work re- 
mains almost constant at zero until oscillations set in at 
around t w 100 (see the inset of Fig. (Ha)). Those oscil- 
lations lead to the slow increase in work source quality 
in the second half of the considered time interval. Al- 
though the oscillations have only small amplitude, their 
work source effect becomes significant due to their fre- 
quency, which is high when compared to the time scale 
of heat dynamics. 

In case (b), W(t, 0) shows strong oscillations of an am- 
plitude orders of magnitude larger than in case (a) from 
the very beginning of the dynamics. Thus, the reason for 
the contradiction to the purity measure result is traced 
back to the problem already touched on in Sec. IVII Al 
Although the purity can be used as a measure for the 
size of incoherent part of the effective dynamics of the 
spin, £° ff , which is associated with the heat flow, it is 
completely insensitive to the size and effects of the co- 
herent par H° s and thus the actual work source effect. 

From this result, we can draw the conclusion that due 
to the fact that the purity measure is only linked to £® 
alone, it can only be used as a necessary condition for 
work source functionality. In order to get the full picture, 
a more detailed analysis of the work and heat fluxes via 
the measure R is necessary. 



C. Work source quality overview 

Computing R for a low photon number parameter win- 
dow (a S [0,4]) and initial spin temperatures ranging 
from to oo (c 6 [0.5, 1]), we find the following trends: 

For a = 0, the overall work source quality of the con- 
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sidered interval t S [0, 200] is generally significantly lower 
than for the corresponding (with respect to c) cases for 
a > 0, and i?(200, 0) ranges from 0.8 to 0.2 with de- 
creasing c. For a > 0, i?(200, 0) takes on values around 
0.9, with a slow increase for higher a and c. The first 
increase can be related to the higher excitation of the os- 
cillator and the resulting bigger amplitude of the position 
expectation. The second trend has to do with a special 
property of the initial states |0)|a) and |l)|a) of which 
the initial state of the xz-SOM is a statistical mixture. 

In order to explain the trend for increasing c, we in- 
voke the first order perturbation theory for the extremal 
initial states |0)|a) and |l)|a) which is applicable to the 
beginning of the dynamics, as long as gt* <C 1 holds with 
g = 10~ 2 /\/2. The calculation is carried out in App. |D] 
Here, we only make use of the result 



U I (t)\Q)\a)=(\0)-iagt\l))\a)+O(g 2 
U I {t)\l)\a) = \ II; 



*y5*|0) I \a) 
igt \0)^L\a)+O(g 2 ) 



(44) 
(45) 



and a* denotes the complex conjugate of a. From this 
form of the state d a * |a), we easily see that in first order 
U 1 (t)\0)\a) factorizes contrary to C/ 7 (i)|l)|a). The purity 
behavior of the initial state p(0) = [c|0)(0|+(l-c)|l)(l|]® 
| a) (a | continuously changes from the |0)|a) case to the 
|l)|a) case. As the |l)|a) state becomes more and more 
mixed into the initial state with decreasing c, a significant 
purity drop happens at earlier times of the evolution. 
The same is true for the heat flow, which is tied to the 
purity drop. With this increased heat flow at the early 
stage of the evolution, Q(200, 0) reaches higher values for 
decreasing c. 

Moreover, the size of work flux for the same change of 
effective splitting of the spin decreases with decreasing c 
until it reaches for c = 0.5. Thus, in the beginning and 
as long as the spin's state is close to its initial occupa- 
tion, the work source effect of the oscillator is reduced or 
suppressed additionally. Clearly, this reduces the work 
source quality and explains the trend seen in the numer- 
ical results. 



VIII. SUMMARY 

Work and heat are related to (thermodynamic) pro- 
cesses, which seem to require external control. In this 
paper, we have argued that work functionality may show 
up in closed bipartite quantum systems, even down to 
the nanoscale. We have shown under what conditions 
the respective subsystem dynamics may be described via 
time-dependent effective Hamiltonians and in this sense 
act as classical driver for each other. 

We have then brought forward the argument that en- 
ergy exchanged under such conditions has to be consid- 
ered as work from the viewpoint of thermodynamics. We 



have applied these results to a simple model confirming 
that a system as small and simple as a single harmonic 
oscillator coupled to a spin can act as a work reservoir for 
the latter. In addition, we have introduced purity based 
and work/heat based work source quality measures and 
discussed their usefulness. We have demonstrated that 
due to the lack of sensitivity to the effects of Hf s , the pu- 
rity based measure is only a necessary condition for work 
source functionality. In general, the implementation of 
a full thermodynamic process within a closed quantum 
system will require driving as well as thermalizing em- 
beddings. 
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Rambach, F. Rempp, J. Teifel, P. Vidal, G. Waldherr, H. 
Weimer and M. Youssef for fruitful discussions. Financial 
support of the Studienstiftung des Deutschen Volkes is 
gratefully acknowledged. 



APPENDIX A: GENERALIZED 
FACTORIZATION APPROXIMATION 

In its original form the FA is formulated for a bipar- 
tite system (cf. [13], Eq. (|3l4p ). Let us now consider a 
tripartite system defined by the Hamiltonian 



H — Hq + Hi + H\2 



Hn 



(Al) 



Let us assume that system has interacted with system 1 
in the past but is now decoupled from system 1. System 
2, however, is supposed to have been uncoupled in the 
past and is now being coupled to system 1 alone. Finally, 
we assume that the combined system 01 and system 2 are 
now in a pure state. We are then left with an initial state 
for the whole system of the form 



p(Q) =poi(0)®/3a(0) 



(A2) 



with p 01 (0) = |*oi (0)) (*oi (0)| and p 2 (0) = 
|* 2 (0))(* 2 (0)|. _ 

We now consider the dynamics of this system with 
respect to the given initial state. In the case that 
tr[p|(i)] ~ 1 holds, the FA is applicable to the whole 
system yielding the two coupled equations 



*fi|*oi (*)> = (#o + #i+ 

+ <*2(*)|Hl2|*2(*)>) |*0l(*)> (A3) 

ih\4> 2 (t)) = [H 2 + (*0l(*)|-ffl2|*0l(*)>) 1*2(0) • (A4) 
By restating Eq. (|A3[) in the form 

ih'm{t)= [^o + -ffi + (*2(t)|Hi2|* 2 (t)),A)i(0] ( A5 ) 

and taking the trace of the Hilbert space of the ancillary 
system 0, we arrive at the result 



ihkit) = Hx + (* 2 {t)\H 12 \$ 2 {t)),pi{t) 



(A6) 
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To get this result, we have made use of the two partial 
trace relations 



tr 
tr 



A 



lr. 



li,B 
)A,B 



= 



A, tr Q B 



(A7) 
(A8) 



Note that in contrast to the case of the FA for a bipartite 
system, the criterion for the applicability of the FA is the 
purity dynamics of system 2 alone. 



APPENDIX B: PURITY OF THE OSCILLATOR 

The purity dynamics of the oscillator in the case of the 
pure a z x interaction can be derived from the solution of 
the Liouvillc-von Neumann equation given in Eq. (|26|) . 



P(t) = 



c\a-(t))(a-(t)\ 



(l-c)\a + (t)){a+(t)\ 
Thus, the oscillator reduced state is 

p„(t) = c|a_(i)> <«-(*) I + (!- c)\a + (t))(a + (t)\ (Bl) 

and taking the square and the trace of this expression, 
we end up with the result for the purity given in Eq. (|27[) . 

P[ P o{t)] = c 2 + (1 - c) 2 + 2c(l - c)\(a_(t)\a+(t))\ 2 . 

For the time-dependent term we find 

\(a-(t)\a + (t))\ 2 = \(a\uHt,0)U+(t,0)\a)\ 2 (B2) 
= \(a\exp(iH_t)exp(-iH + t)\a)\ 2 (B3) 

with H± = H ±Ax±Hf-. Making use of the displacement 
operator D{a) = exp(aa^ — a*a) and its properties 



D(-a)xD(a) 

D(-a)pD(a) 
we can express H± as 



muj 



-Re(a) 



p + V2mw Im(a), 



H± = D(-l3 ± )H D(p ± ) 
and therefore have 



C 



U±(t,0) = e- l0t D(-[3±)eM-iH t)D(P±) 



(B4) 
(B5) 

(B6) 
(B7) 



with (3± = qpA/ ^2muj3 U p to constant factors or a phase, 
respectively, which are irrelevant for the computation of 
the modulus in Eq. (|B2|) . With the help of the relations 

b(a)b{P) = exp[ilm(a/3*)]b(a + 0) (B8) 
cxp(— iH t)\a) = exp(— iuj t/2)\aexp(— iu t)) (B9) 
we arrive at 

\a±(t)) =exp[i(f>±(t)]\{a+/3±)exp(-iu t)-l3±}. (BIO) 
Finally making use of the relation |(a|a')| 2 = exp(— \a — 



a'\ ) yields the result 
\(a-(t)\a+(t))\ 2 =exp 



■ sin 



-U! t 



(Bll) 



APPENDIX C: AMPLITUDE OF THE SPIN'S 
EFFECTIVE SPLITTING IN THE 2-SOM 

In the case of the z-SOM and applicable FA, the effec- 
tive Hamiltonians of the spin and oscillator are found to 
be [Eqs. (150121)1 ] 

Hf(t) = (± + \{x)(t))a z (CI) 



Hf = [Ho + A(l - 2c)xj (C2) 
and the latter may be rewritten in the form 

Hf = hjj a (x 2 + P 2 ) + const. (C3) 



with the dimensionless position and momentum opera- 
tors 



X 



P = P, 




(l-2c) 



(C4) 
(C5) 



where X — ^/muj x, P = p/^/mu) . Making use of 
the properties of the displacement operator b{a) in 
Eqs. (|B4IB5|) . one finds that 



Hf = D(- 7 )if £>( 7 ) 



with 



7 



/ A 2 
2muu, 



(l-2c). 



(C6) 



(C7) 



In order to compute the peak-to-peak amplitude of the 
effective spin splitting 

A^f = A((£) max -(£) min ) 

we need to evaluate 



(C8) 
(C9) 



(x)(t) = (a(t)\x\a(t)) 

= (a\exp(iHft)&exp(-iHft)\a) 

= (a|£>(-7) exp(iH t)D(j)xb(-j) exp(-iH t)D(j)\a) 

= (a + i\exp(iH i)xexp{-iH t)\a + i) +7, (CIO) 

where we have used Eqs. (IB4IB8IC6P assuming ael and 
defining 7 — y / 2/(mu )'-f. Now, we can see that the first 
term is just the time evolution of the expectation value of 
the position of the original oscillator described by H Q for 
a coherent initial state |a+ 7 ). With the help of Eq. jB9j 
it is straightforward to show that (A~) max — (A) min = 
2V2|q + 7 1 and therefore 



Awf = 2A 



-|a + 7|. 



(Cll) 



Note that this result is only exact if the Hamiltonian gov- 
erning the oscillator's dynamics is Hf and Cf \p(t)] = 0, 
that is if the FA is exact. Still, if the FA holds in good 
approximation, Eq. (jClip is a good approximation as 
well. 
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APPENDIX D: FIRST ORDER 
TIME-DEPENDENT PERTURBATION THEORY 
FOR PURE INITIAL STATES OF JCM 

It is convenient to apply the perturbation theory in 
the interaction picture. All interaction picture quantities 
are denoted by a superscript "7" . The expansion of the 
time-evolution of the state is given by 



\* I (t)) = \*(0))+U[(t)\*(0))+O(g 2 



(Dl) 



and the first order contribution to the time-evolution op- 
erator U((t) is given by (see, e.g., [HI, p. 207ff) 



U[(t) 



i J AtV 1 ^) (D2) 



and 

= gexp[i(H s + H )t](& + a + a— a') exp[— i(H s 



H )t] 
(D3) 



is the interaction operator in the interaction picture. Ac- 
cording to the RWA, only terms of the interaction are 
kept which are time-independent in the interaction pic- 
ture, thus the last equality. From Eqs. (|D3[) and (|D2|) . it 
follows that 



U((t) = -igt(a+a + <7_a t ). 



The time evolution of a state in the JCM is therefore 
given in first order perturbation by 

\* z {t)) = [l-igt(&+& + |*(0)) + 0(g 2 ). (D5) 



Here, we consider |0)|a) and |l)|a) as initial states. 
Together with 



a* la) = 



d a* 
da~* + T 



(see [23]) we find for those states 



U I (t)\0)\a)=(\0)-iagt\l))\a) + O(g 2 ) 
U 1 (t)\l)\a) = (\l)-i^gt\0)^ \a) 

-igt\0)-^\<x) + O(g 2 ) 



where 



da" 



d 

da* 



exp 



2 



E 

71=0 



a" 



(D4) (a* is the complex conjugate of a). 



(D6) 



(D7) 
(D8) 



(D9) 
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